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Abstract: We investigate the effect of adding a Chern-Simons term coupled to an 
axion field to SU(2) Einstein- Yang-Mills in a fixed AdS^/Schwarzschild background. 
We show that, when the axion has no potential, there is a phase transition between 
a Reissner-Nordstrom black-hole and one with a non-abelian condensate as per the 
vanishing Chern-Simons case. Furthermore, by giving the axion field a mass, one 
observes a phase transition between a Reissner-Nordstrom black-hole with axion 
hair to a "superconducting" phase which also has a non-trivial axion profile. We are 
able to perform a preliminary analysis for this interesting case and observe that we 
can shift the critical temperature at which the phase transition occurs and observe 
interesting features of the order parameter scaling form. 



1. Introduction 



In recent years the subject of phase transitions in black holes has been extensively 
studied. Much is to be attributed to the development of AdS/CFT [1] within 
the context of condensed matter systems, especially in holographic superconductiv- 
ity [2], [3] (see [4] [5] [6] for comprehensive reviews on the subject). In such models 
superconductivity is described by a phase transition of an asymptotically AdS black 
hole which admits both a Reissner- Nordstrom (RN) and an ArfS'/Schwarschild solu- 
tion. In the initial proposal an abelian condensate is described by a charged scalar 
field acquiring a VEV on the boundary of the AdS space or, in the dual gravita- 
tional picture, the black-hole developing scalar hair. The model was extended to in- 
clude phase-transitions of AdS black-holes involving non-abelian condensates [8] and 
these led to phenomenologically promising models of p-wave superconductivity [9]. 
Since then, models of p-wave holographic superconductivity have been widely stud- 
ied [11] [12] [16] [17] [18] [19] [20]. 

In [21] it was shown that one can also observe characteristics of Chern-Simons 
(CS) interactions in superconductivity described by abelian condensates by coupling 
a CS term to a neutral axion field in the four-dimensional Einstein-Maxwell action. 
In this case the condensate is still described by an external scalar field which doesn't 
couple to the axion field, hence the condensate profiles are identical to those with 
a vanishing CS term. Vortex solutions of the system lead to properties of pure CS 
systems such as the magnetic field peaking outside of the core of the vortex. In this 
paper we wish to make progress towards including CS effects in four dimensional AdS 
black-hole phase transitions involving non-abelian condensates. Our mechanism will 
be very similar to the previously mentioned case: we will couple an £77(2) CS term 
to an axion field, solve the system in the bulk and project it to the boundary where 
we hope to observe interesting features of the dual field theory. In this case however, 
since the condensate is provided by the gauge field itself, we are effectively coupling 
the axion field directly to the condensate through the CS term and this will lead to 
interesting novel observations on the profile of the order parameter and the space 
of possible solutions describing the phase transition. Unlike the case of the abelian 
condensate we don't expect to observe CS characteristics in the boundary field the- 
ory itself, our approach here is to investigate the holographic effects on this theory 
observed via the CS coupling in the bulk. The hope is that through this study one 
could in the future investigate the full effects of the CS term on a model of p-wave 
superconductivity. 

The paper is organised as follows: in section 2 we introduce the system consist- 
ing of £77(2) Einstein- Yang-Mills with the inclusion of an axion field coupled to 
a CS term. We make an ansatz for the gauge field, derive the equations of motion 
of the system and by expanding the fields at the AdS boundary discuss the relevant 
dual thermodynamical variables of interest to the problem. Section 3 is devoted to 
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numerical solutions to the above system. There are two parts: the first in which the 
axion potential is set to zero and the second in which we switch on a mass term for 
the axion field. We show here that these lead to very different results. In section 4 we 
calculate and evaluate the free energy of the different phases by using the AdS/CFT 
dictionary and in section 5 we provide short analytical results of the T ~ T c region 
which give useful general behaviours of quantities of interest. Finally in section 6, 
we provide a brief summary of the results obtained and point the reader to directions 
for future work. 



2. The System 

Our starting point will be SU (2) Einstein- Yang-Mills with a Chern-Simons term 
coupled to a dimensionless axionic field 9 in the AdS±/ Schwarzs child (AdS/Sch) 
background. 



S= / dx A y/=G[R-A+\c) (2.1) 



2k\ J V V 

where 

C = --Tr (F^F^) + -^L=9Tr {F A F) + \ 2 d tl 9d»9 + V{9) (2.2) 
4 y—G 

where F° v = d^A^ - d v A° + e^A^A^, the cosmological constant A = -j^, L is 
the AdS radius, k\ = 8itGn, A is a dimension one constant (usually called the axion 
decay constant) which we set to 1 manually and V(9) is the axion potential which for 
the moment we will leave unspecified, k is a constant useful in keeping track of the 
relative Chern-Simons contribution to the action. It is important here to note that 
this set-up manifestly violates parity symmetry under the transformation 9 — > —9. 

We work in the non-backreacting limit of large g where we can take the AdS/ Sch 
ansatz for the metric 



ds 2 = ^ 



- 1 - ^ dt ' + dx ' + dy z 



L dr 2 

+ ^IT4' (2 ' 3) 



in which indicates the position of the black hole horizon. We take the gauge 
field A = A a ^r a dx^ where r a are the generators of the SU(2) algebra such that 
[r a , r b ] = e ahc r c . The effective boundary Chern-Simons coupling term is k9 evalu- 
ated on the boundary and we must impose that this is quantized to work with a 
reasonable Chern-Simons theory (see [7]), this will put restrictions on the values of 
9 at the horizon as will be shown later. 
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In [8] it was shown that, for the case of k = a p + ip gauge field ansatz 

A = (f)T 3 dt + u {r l dx + r 2 dy) (2.4) 

leads to a second order phase transition from a RN anti-de Sitter black hole 
to one with non-abelian condensates, with non zero lo acting as the condensate. 
However in [9] the authors argue that these backgrounds are unstable to small per- 
turbations that turn them into the less isotropic p-wave backgrounds, at least close 
to T = T c . It is an important question whether with the introduction of a dynamical 
axion term the p-wave backgrounds remains thermodynamically preferred, however 
we will not study this here and focus on the interesting effects of the axion on the 
p + ip transition with the hope to report on the above problem in future work. 

Throughout this paper we decide to work in units of the AdS radius, L = 1 and 
the horizon radius — 1. Hence in this case £, ^ — > r is a dimensionless parameter. 
Similarly, we can form dimensionless fields by coL — > Co and 4>L — > 0. Hence bear 
in mind that u>, and r are all dimensionful, as are all other fields in the theory, 
however the results presented in the paper are for dimensionless combinations of Co, <p 
and — with r^, L — 1 manually. 

The equations of motion one derives from the action are 

f + V - 4n 2 M 0^ 2 + - 2 k9'co 2 = (2.5) 
r r 4 (l — ^) r 2 

r ( r 3_i) r 4 (l-^) r\l-±f V r 2 (l-^) V 1 ; 

if - '^V-V IP® - *L- { *S)' = (2.7) 
r(r 3 -l) r 2r 2 (l - ^) r*(l - ±) W } 1 ; 

where in the above ' denotes differentiation w.r.t r. Note that as per [21], k6' 
acts as an effective Chern-Simons coupling in the boundary theory. The fields have 
asymptotic behaviour (again these should be thought as dimensionless fields at the 
horizon with r replaced by ^ everywhere) 



= <Mr-l) + 2 (r-l) 2 + ... (2.8) 
oo = oj + u>i(r - l) 2 + ... 
= *o + *i(r-l) + .... 

at the horizon and 



4> = Po + - + ... (2.9) 

r 
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Wi 

u = — + .... 

r 

6 = 6 + ^ + .... 

at the asymptotic boundary at large r, where the expansion for 9 is at m = 0. 
Through the gauge/gravity correspondence we can associate thermodynamic quan- 
tities to the above variables (2.9) in the expansions of the fields at the asymptotic 
boundary. The holographic dictionary states that a bulk field is dual to a field theory 
operator on the boundary through the relation 

(exp J O O) = exp{-S os [<f> ]), (2.10) 

where S os denotes the on-shell action. More precisely, the boundary value of 
the bulk field <fi acts as a source for the corresponding field theory operator on the 
boundary. The method of calculation of the precise relationships between thermody- 
namic quantities of interest and bulk fields is illustrated fully in Appendix A of [10], 
to which we adhere. Therefore the charge density 

p = w/ drAtLy = W9 Pu (2 ' n) 

and the chemical potential 



V = WJ^zMbry = T^lZPO- (2.12) 



— A I — 
2kjg tlbry ~ Iklg 1 

note that here the asymptotic field values are untilded, hence they are dimen- 
sionful. In order to remove the dependence on /c 4 from these quantities we make the 
redefinition p = j^sp where the extra factors of tt are included for numerical con- 
venience. Similarly the AdS/CFT correspondence maps the Gibbons-Hawking [22] 
temperature Tqh = ^ of the black hole to the Temperature T of the CFT and, as 
explained in detail in [5], for a scale invariant theory at finite temperature in equi- 
librium there is no other meaningful scale which we can compare the temperature 
to, i.e. all non-zero temperatures should be equivalent. Hence, it is only physically 
meaningful to compare scale invariant quantities formed from the bulk fields pro- 
jected on the boundary, to this purpose we define the effective temperature of the 
dual CFT as which from the above definitions is easily seen to be equal to 



As mentioned below equation (2.4), a non- vanishing value of to at the boundary 
acts as a condensate for the dual field theory. Hence 



2k 2 g ~r-^ bry 2k 2 g VV l 2g 



J = -^2-J 2 drA x \ bry = ^j-m = -Wi (2.14) 
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will be our order parameter, where in the last line we have made a similar 
redefiniton of u to remove the dependence on k\. As per the effective temperature 
defined previously we construct a physical scale invariant order parameter i, which 
evaluates to 

i - -A (2.15) 
P Pi 

In similar fashion, the bulk axion 9 maps to a pseudo-scalar field theory operator 
on the boundary field theory, with 9q acting as its source. The boundary behaviour 
of the axion field is such that 

9(r) = 9 r- A - +9 ir ~ A+ (2.16) 

where A± indicate the dimensions of the corresponding boundary operators in 
the dual field theory. Note that as per a scalar field in AdS space, the dimension of 
the boundary operator is dependent on m, the mass of the corresponding field in the 
bulk. In the case of m = then the asymptotics are those shown in 2.9, however for 
the case where m^fl then the analysis of varying m shown below, which is for the 
case of 9q not held fixed in the numerical procedure, corresponds in the field theory 
perspective as variations of the dual operator dimensions. 

To determine solutions to the equations of motion (2.5) we adopt a numerical 
"shooting" procedure, which involves allowing for a non-vanishing constant term 
Wo in the expansion for Cj at the boundary and then manually ensuring that this 
vanishes 1 by carefully changing the choices of Cjq and <j>i. For the case of vanishing 
axion these three conditions are enough to determine a one parameter family of 
solutions which seeds the values at the asymptotic boundary, as explained in [9]. In 
the case of a non-vanishing axion we have three second order differential equations 
and thus six integration constants. The constraints are normalizability of u, the 
values of 4>i and uio an d the values for to and t\, hence we have six constants and five 
constraints which makes a one-parameter family of solutions as per the previous case. 
We proceed in the numerical analysis working at constant po (the grand canonical 
ensemble) and enforce that the asymptotic value of the axion 9o is constant manually. 
For the case of a non-vanishing axion potential in which the asymptotics for 9 cannot 
be controlled manually we proceed under the assumption that varying m smoothly 
can enforce this condition, it is however not what has been done here and therefore 
this section should be taken only as a preliminary investigation on the effects the 
potential may have on forms of solutions with w,k^0. 



^^This is the normalisability condition for ui. 
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3. Solutions 



This section is devoted to analysing numerical solutions of the above system. There 
are two variables which can be tuned by hand : V{9) and k both subject to the overall 
solution being normalizable and thermodinamically preferred (see section below). We 
will start with the analysis for a vanishing axionic potential and proceed, in the next 
section, to include a mass term for the axion field. 

3.1 V(0) = 

For the simplifying case of V(8) = we will observe the effect of raising the Chern- 
Simons parameter k on the space of possible solutions of the system, the profile for 
the axion field and the form of the condensate. The normalizable perturbations to 
the gauge field u have many solutions with increasing nodes. We will restrict to so- 
lutions of the form shown in Figure 1 where a) is a monotonically decreasing function 
with no nodes as these are believed to be thermodynamically favoured over the other. 




Figure 1: Thermodynamically preferred form of the normalizable perturbation uj in which 
this is a monotonically decreasing function with no nodes. 

Throughout this section we work with to = 1 and t± = 0, we also drop the ex- 
plicit tildes from the field values, bear in mind however that the quoted results are 
for dimensionless tilded forms of the fields. Figure 2 illustrates all possible values (to 
within the numerical accuracy of the procedure) for ujq and 0i which yield normal- 
isable solutions with a non-vanishing condensate. The normal phase corresponds to 
ojq = whilst the phase with non-zero uq will be referred to as the superconducting 
phase. The blue/leftmost line corresponds to the case of k — 0, where the above 
ansatz for the gauge-field was shown to cause a second order phase transition be- 
tween a RN black-hole and one with non- vanishing non-abelian condensate. Curves 
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to the right of this are for increasing values of k in steps of 0.04. A maximum value 
of k — 0.2 is reached before the numerical procedure breaks down, at this value the 
CS term is large enough that one cannot ignore the back-reaction of the gauge fields. 
Note that in this scheme where V(6) = we see that all curves tend to the same 
value of 0i as u — > 0. In this case one finds that pi is never greater than about 3.71, 
which leads to a constant (k independent) value of the critical temperature 

T c ^\mJgL P . (3.1) 

We will show in a later section that the independence of T c on k is justified 
analytically. 



oj 




Figure 2: The parameter space of possible solutions to the non-linear problem. The 
normal phase is described by oj = whereas points on the line with non-zero uj describe 
the superconducting phase. As charge density and chemical potential are held constant, 
Temperature increases towards the left. The curve to the left is k = 0, then increasing k 
to the right in steps of 0.04. Hence the highest k shown in black is k = 0.2. 

As is seen in Figure 2, for a narrow range of k the transition switches to first 
order. This is best seen in Figure 3 where the same plot is presented for k = 0.12, 
showing that the transition is first order, i.e. as one decreases the temperature for 
non-zero u there is a non-continuous phase transition with a corresponding jump in 
the free energy. Increasing k even further restores the transition to second order. The 
first order phase transition is seen for values of k ranging between k = 0.115 to 0.123. 

The axion profile is shown in Figure 4. For apart from a sharp rise around 

the position of the horizon we see that the axion is a constant everywhere. The 



-7 - 




3.700 3.705 3.710 3.715 



Figure 3: First order phase transition at finite k = 0.12. The "bump" seen extending 
beyond the region where the two phases meet at oj = signifies that as one lowers the 
temperature there is a finite jump in the free energy between the two phases, and hence a 
first order phase transition. 

different curves correspond to increasing values of k, with the blue line corresponding 
to k — 0. The overall shape of the axion seems independent of k. The fact that 9 
is a constant on the boundary does not mean that the black-hole develops axion- 
hair. One can effectively remove this by simply adding a constant term to the F AF 
term in the action. Solutions with the same asymptotic 9q values for the axion field 
were compared here. Since in the case where V{9) = a constant axion profile is a 
solution to the equations of motion, the axion contribution disappears entirely from 
the normal phase, which is simply AdS — RN . In this case we are truly comparing 
a phase with no axion to that with an axion, this is not the case when we switch on 
a potential term as is discussed below. 

Finally, in Figure 5 we show the effect of different values of k on the condensate 
at the boundary. The plot shows the condensate as a function of |r, with the highest 
line corresponding to the K = case, then increasingly lower results at increasing k 
in steps of 0.04. The condensate suffers from a suppression as k increases, the form 
of which is investigated analytically in a later section. 

3.2 V{9) ^ 

In this section we switch on a mass for the axion field so that V(9) = m 2 \ 2 9 2 — > m 2 9 2 
and we will investigate what happens to the phase transition as one varies the mass 
m at finite non-zero small n, we restrict the analysis to k — 0.08, below the region 
where we see a first order phase transition. Note that with the contribution from a 
potential term a constant axion profile is not a solution to the equations of motion 
(2.5), as is indeed seen in Figure 7 even for the simple case of k — 0. This implies 
that we have a non-zero contribution from the axion kinetic and potential terms 
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Figure 5: The condensate profile for increasing k. The highest line corresponds to k = 0, 
then increasing k downwards in steps of 0.04. There is a small suppression in the condensate 
profile for increasing k. 

present in the "normal" phase where u — 0. The comparison is thus between an 
AdS'-Schwarzschild black hole with a non-trivial profile for the axion to a "supercon- 
ducting" phase with a condensate and an axion. Whilst the condition of comparison 
of the same asymptotics for the 9 field is easily met in the case of vanishing poten- 
tial, this cannot be done straightforwardly here. In this case the asymptotic value of 
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the axion field varies as I vary the mass m, this is because the profile of the axion 
varies as a function of m and therefore the extracted values of 8q will be different for 
each profile of the axion in the case where this is not held constant in the numerical 
procedure, which applies here. Hence comparison with the same asymptotic for both 
phases requires smooth variations of m with u> so as to maintain 9 constant. This has 
not been done in this section, which rather illustrates how varying m values changes 
the form of solutions with differing values of 9 with constant u and constant k. It 
remains an important question which solution is thermodynamically preferred when 
the axion has a potential term contribution, whether a phase transition occurs in 
the first place and whether this term has noticeable effects on the space of solutions 
of the system. These are the issues which this preliminary investigation wishes to 
address. 

In Figure 6 we show the effect of raising m on the space of possible solutions 
of the system. The blue/leftmost curve corresponds to the case of vanishing m = 
and curves to the right of this are for increasing values of m in steps of 0.1, the 
final curve is m = 0.45 as already for m = 0.5 the back-reaction cannot be ignored 
and the numerical procedure breaks down. We see that increasing the mass of the 
axion has the effect of shifting the phase transition curve to higher values of 0i as 
w — > whilst preserving its shape. This has the interesting effect of lowering the 
critical temperature T c at which the transition takes place. Unfortunately we are 
restricted from investigating the region of large m from the numerical procedure. 
With variations in m we also observe a variation in the shape of the axion profile. 
This is shown in Figure 7 where, contrary to changing k, the axion has a non-trivial 
profile in the bulk and changing m doesn't correspond to a simple shift for 6. 

It is interesting to observe the effects of m on the order parameter of the field 
theory. We saw in the previous section that when m = 0, varying k had (to within 
the tested numerical range of parameters) a small effect on the T ~ T c region of the 
order parameter. In Figure 8 we have plotted the order parameter against ^ for 

J- c 

curves with different values of m. The blue/highest curve is the m = case and the 
curves below this are for increasing m. Each curve is plotted with its corresponding 
value of T c . These results are all for the choice t = 1. Given that t appears coupled 
to m as the Chern-Simons interaction term on the boundary it is evident that fixing 
m^O means that varying t has analogous effects to the system to varying m. 

4. Free Energy 

In this section we make use of the AdS/CFT correspondence to calculate the free 
energy density for the superconducting phase /. We are interested in the scale in- 
variant quantity 44 — ^~{% N , the difference in free energy densities between the 
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Figure 6: The effect of the mass m on the space of possible solutions at k = 0.08. The 
blue/leftmost curve corresponding to m = then increasing m towards the right in steps 
of 0.1, the final curve is m = 0.45. 
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Figure 7: The effect of the mass m for the axion field at k = 0.08. Increasing m upwards 
in steps of 0.1, the blue/lowest curve corresponds to m = 0. 

normal and the superconducting phase. If this remains negative, then the symmetry 
breaking phase is preferred. 

By the AdS/CFT dictionary the free energy is given by the on-shell Euclidean 
action with appropriate counter terms to cure divergences. Therefore we work with 

2k\S = S grav + SMaxw + ScS (4-1) 

where 
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Figure 8: Increasing m effects on condensate at fixed k = 0.08. The highest/blue curve 
corresponds to m = 0, the lower curves represent increasing m in steps of 0.1. There is a 
significant effect on the profile of the condensate. 



S g rav = I d A x^ (R + 6)+ I d 3 x^ (-2K + 4) , (4.2) 

*J J foe 

the second term is the usual Gibbons-Hawking [13] boundary term needed to 
have a sensible variational principle, 7 is the induced metric on the boundary and 
we also add a boundary cosmological constant term to regulate the action, then 

S Maxw = ~ j d 4 x^Tr (F^Fn (4.3) 

is the usual Maxwell term with no further counter-terms needed as we assume 
that the gauge field goes to zero at the boundary sufficiently quickly, finally 

S cs = J d A x^j {dpdm + V{6)) + J d 4 xn9Tr (F A F) (4.4) 

where these terms are the usual kinetic, potential and Chern-Simons terms for the 
axion field. We then proceed to evaluate this action on-shell with the Euclideanised 
t — > ir metric to obtain 



1 —( ( f ) uj) 2 + 8iK9(cu 2 ( l>y + r 4 (l- , -^)(e') 2 + r 2 V(9) (4.5) 



2r 2 ( 1 - ^ > V r 



where the first term is given by the normal phase Maxwell term, f^ s denotes 
the CS terms evaluated in the normal phase and the remaining terms are obtained 
from the numerical solutions. The terms in f^ s are simply the kinetic and potential 
axion terms in the normal phase, whose presence or disappearance is determined 
non-trivially as they don't have direct couplings to u (which vanishes in this phase). 



From the equations of motion (2.5) we can see that in the case of vanishing axion mass 
m = then the kinetic term in the normal phase, which is proportional to 9' doesn't 
contribute as a consistent solution for the axion is simply to have a constant profile. 
Hence in this case, f® s = and we don't have to worry about CS contributions 
to the free energy coming from the normal phase. We now proceed to evaluate the 
remaining terms in order to obtain the overall change in free energy. Let's start with 
the term involving k, this is 



and hence when we Euclideanise the action, this term acquires a factor of i from 
the dt component of the gauge field. This means that the term is irrelevant in the 
partition function e~ s and thus doesn't contribute to the free energy. This remains 
true in the case of a non-zero axion potential. 

The kinetic term for the axion field contributes to the free energy density in the 
form 



which contributes a finite amount to the free energy, changing with k as the 
corresponding solution for 9 changes. It might at first sight appear surprising that 
this term remains finite over the integration due to the direct coupling to a factor of 
r 4 which diverges at infinity, however we can see from the asymptotic behaviour of 
the axion 2.9 (9') 2 m r~ 8 that this factor is cancelled in the overall integral, which 
thus remains finite. A sketch of the profile of the derivative of the axion term is 
shown in Figure 9, note how the axion decays to zero for large r. Figure 10 shows 
the effects that increasing k has on the free energy through the changes in profile 
of the axion term. For the values of k tested in the numerical procedure, the curve 
remains negative and thus the axion phase is preferred. 



In the case of non-vanishing potential things become greatly more complicated. CS 
contributions are present in the normal phase as 9 having a constant profile is no 



4.1 V{9) = 




(4.6) 




(4.7) 



4.2 V{9) ^ 
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Figure 9: The numerical solution for & . At large r the profile is vanishing, this illustrates 
how the kinetic contribution of the axion term remains well behaved despite having explicit 
couplings to positive powers of r. 




Figure 10: The effect of increasing k on the free energy density. The lowest blue curve 
corresponds to k = then increasing k upwards in steps of 0.04. For all k, Af remains 
negative. 



longer a solution to the equations of motion. Hence in this case the normal phase 
contributes 

f No = -J + f dr (l - ^) i ') 2 + r2y (0?) U (4-8) 

and we have to make sure we consider the resulting kinetic contribution arising 
between the difference in kinetic axion terms in both the normal and superconducting 
phases. Most importantly the term 

/•oo 

fry) = / drr 2 V(6) (4.9) 

is divergent as V(0) = m 2 6 2 and 9 is non-vanishing at asymptotic r. One can 
hope that the divergent contribution from the normal phase here can cancel that of 
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the superconducting phase, however this cannot be the case exactly since the two 
phases have different profiles for the axion 9. Therefore to make progress we could 
proceed in two distinct ways: by a systematic holographic renormalization procedure 
(see [14]) in which the necessary counterterms are computed manually and the diver- 
gences are cancelled explicitly or by ensuring that the numerical procedure is such 
as to maintain a constant value for # a t the boundary, in which case the divergences 
should cancel exactly in both phases (clearly both methods should eventually agree). 
This section is intended as an introduction to the case of V(8) 7^ and thus we won't 
carry out the manual renormalization procedure here. Therefore, it remains an open 
and important questions whether a non-zero mass m for the axion can de-stabilise 
the "superconducting" phase or not. We wish to report on this (and especially on 
the second numerical method for cancelling the divergences) in future work. 

5. Analytic Calculations 

Even though the set of equations (2.5) require numerical solutions, one can still 
obtain an analytical feel for the behaviour of thermodynamic quantities of interest 
by looking at the region close the phase transition. In this region the fields are small 
so one can trust the series solutions for the fields obtained by matching the series 
expansions at the horizon and at the boundary, the procedure is detailed in [15]. Even 
though it is an excellent method to gain analytic understanding for the behaviour 
of quantities of interest this procedure is only an approximation valid in the region 
of small fields, hence one should be sceptical about the precise numerical results 
presented. By inserting equations (2.8) into (2.5) and matching to the expansions at 
the boundary we obtain for the case of m = that 



where C is an arbitrary proportionality constant chosen in the analytical proce- 
dure. By choosing C = 1.95 we can recover the numerical result 3.1. So the critical 
temperature doesn't depend on k as is observed in Figure 2, and 



which shows the dependence of the condensate on k. The precise numbers ap- 
pearing in this relations are unimportant, it is intended to demonstrate the explicit 
behaviour with varying k. It also indicates that the transition at finite k has simple 
critical exponents given that T c is k independent. This should be trusted close to the 
transition in the vicinity of T = T c . From Figure 5 we can see that this analytical 
behaviour matches well with the numerical results. Switching on a mass term for the 



T c = CyfgLp, 




J 



P 




(5.2) 
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axion gives a complicated expression for T c (m) which we won't include for simplicity, 
however it is important to see that even analytically when we include a non vanish- 
ing potential for the axion field then the critical temperature shifts as a function 
of the mass, as observed in Figure 6. Furthermore, this indicates that in this case 
the transition may not have the simple ^Jl — ^- behaviour close to the transitions, 
depending on the value of m. 

We must also impose that the boundary Chern-Simons term is quantized. The 
requirement 

Ti 

K0 r „ = - (5.3) 

is translated to a requirement on the series expansion at the horizon for 9 by the 
analytical procedure. We obtain that 

n 96tn 
to = 16^ " 7 + 192^' (5 ' 4) 

this makes sense in this approximation where T is close to T c and t ~ 0. Working 
with to — 1 in the numerical procedure amounts to a rescaling of the analytical to 
by 16-kk and working at n — 1. 



6. Conclusions 



In this paper we have investigated the effects of adding a Chern-Simons term, cou- 
pled to an axion field, to the phase transition between AdS Reissner-Nordstrom and 
AdS black-holes with a non-abelian condensate. In particular we considered the two 
cases of vanishing potential for the axion field and providing this with a finite mass 
V(9) = m 2 X 2 2 . In the first case we showed that increasing k led to a transition 
which alternates between first order and second order. As oo — > all solutions ended 
at equal 0i and no change in T c was observed. In terms of the order parameter 
we showed that this was suppressed with increasing k, even though for the range 
of k tested the effects were small. We showed analytically that the condensate has 
a (a + bn 2 )~2 for constant a,b dependence which matches well with our numerical 
results. There were no noticeable changes in the results when varying the choice for 
the expansion of 9 at the horizon. With calculations of the free energy we showed 
that the superconducting phase remains dominant below T c . 

When the axion field is given a finite mass things are significantly different. 
Firstly, the axion field develops a non-constant profile in the bulk which makes the 
asymptotic comparison of the two phases difficult, the space of solutions which admit 
a non- vanishing uj is shifted to higher values of 4>i as oj — > 0. This in turn has the ef- 
fect of lowering the critical temperature T c at which the transition takes place. Being 
non-constant the axion contributes to the normal phase as well, hence in this case 
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the transition is between an AdS/ RN black hole with some axion hair to a "super- 
conducting" phase with axion hair as well. Also, as pointed out inside the paper, the 
asymptotic matching of the axion is hard when the potential is present and therefore 
the results shown are only for constant u and k but varying 9 . Indeed it is crucial 
to extend the algorithm used for the results of this paper to one which can maintain 
constant 6 values at the boundary. Even though this might appear a straight-forward 
task conceptually (after all it is the addition of one further constraint to the numeri- 
cal solver) the implementation in terms of the numerical procedure presents problems 
due to the highly non-linear nature of the axion interaction terms. To make a true 
comparison between phases in which u = to those called "superconducting" one 
would have to maintain the same boundary asymptotics for 8 which requires smooth 
variations of m as u is varied. This is an important task which certainly should 
be investigated further. Nevertheless, the presence of the potential has significant 
effects on the characteristics of the phase transition at constant u, suppressing the 
order parameter. Similarly, one would like to show that for all regions of m the free 
energy contribution of this potential term is small enough so as to not make the 
overall energy change positive, signifying that the transition is not preferred. Once 
the constant asymptotic values for the axion are maintained the divergences should 
cancel in both "superconducting" phases and the results should be reliable. We wish 
to report on this area of work in future work. 

In conclusion we have observed that a CS term has important effects on the phase 
transition of AdS black-holes with non-abelian condensates. It would be desirable 
to investigate the region of large m and large k, where in the most optimistic case 
(with high-precision numerics) one could also push T c down to zero. As is apparent 
from the paper, the most important obstacle in making progress with this system is 
the complexity of the numerical procedure. One would really like to be able to probe 
the limiting regions of the parameter space. For example we would like to investigate 
the interesting region of large k where the system is indicative of admitting a region 
of non- vanishing u for large <f>i. The numerical procedure is reliable up to the point 
when terms become large enough to back-react significantly on the geometry and the 
approximation breaks down. 

One very important direction for further research involves modifications of the 
gauge field ansatz. The next step is to investigate the effects of the CS term in 
the context of p-wave Holographic Superconductivity [9], this requires a less istropic 
form for the gauge field. One way to make progress is to include constant values in 
both directions of the gauge field ansatz, i.e. A ps u(ar } dx + 6t 5 l dy) for a, b constants 
and vary these numerically. Then one can hope that in regions where a » b we 
could drive the system to resemble one which has fewer isotropic characteristics. 

Regarding the above system, one immediate desire is to push the analysis to the 
fully back-reacting case of finite g in which one would observe the complete effect of 
the CS term. Furthermore, one can freely change the profile for the axion potential 
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V(8), for example one could include a X8 A term, this may lead to novel interesting 
features not observed here. In general, the axion potential system deserves detailed 
analysis with regards to the free energy and the asymptotics, where all scaling sym- 
metries of the equations of motion could determine whether a comparison between 
the two phases is reliable even without explicit matching of 9 . The OF A F term 
can arise, for example, in anomaly cancellation in string-theory. It remains an open 
and interesting problem to find a suitable string-theory reduction which leads to the 
above set-up. 
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